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Can the early X-ray afterglow of GRBs be explained 
by a contribution from the reverse shock? 
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ABSTRACT 

We propose to explain the recent observations of GRB early X-ray afterglows 
with SWIFT by the dissipation of energy in the reverse shock which crosses the 
ejecta as it is decelerated by the burst environment. We compute the evolution 
of the dissipated power and discuss the possibility that a fraction of it can 
be radiated in the X-ray range. We show that this reverse shock contribution 
behaves in a way very similar to the observed X-ray afterglows if the following 
two conditions are satisfied: (i) the Lorentz factor of the material which is 
ejected during the late stages of source activity decreases to small values 
r < 10 and (m) a large part of the shock dissipated energy is transferred to a 
small fraction (£ < 10~ 2 ) of the electron population. We also discuss how our 
results may help to solve some puzzling problems raised by multiwavelength 
early afterglow observations such as the presence of chromatic breaks. 

Key words: gamma ray: bursts; shock waves; radiation mechanisms: non- 
thermal 



1 INTRODUCTION 

The X-Ray Telescope (XRT) on board the SWIFT satellite has for the first time allowed a 
follow-up of the X-ray afterglows of GRBs starting within one minute of the BAT trigger 
(Burrows et al. 2005a). These early afterglow observations have revealed several surprising 
features which cannot be easily understood in terms of the usual interpretation where the 
afterglow comes from dissipation in a forward shock propagating through the source envi- 
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ronment. At very early times immediately after the burst prompt emission, the afterglow 
first exhibits a steep decrease of temporal slope a± ~ 3 - 5 (Fx oc t~ a ) (Tagliaferri et al. 
2005). It is often followed by a much shallower part with 0.2 < a 2 < 0.8 which can last 
for several hours until a more standard slope 1 < a 3 < 1.5 is finally observed (Nousek et 
al. 2005). Moreover flares with sharp rise and decay times are often present, superimposed 
on the power- law evolution (Burrows et al. 2005b). In most cases, the spectrum remains 
essentially constant through the breaks which may indicate that a single physical process is 
responsible for the whole X-ray emission. The forward shock could be such a process but it 
seems unable, at least in its simplest version, to account for the early slopes a\ and a 2 - It 
has been suggested that the shallow part of the light curve could still be produced by the 
forward shock if it is continuously feeded in energy by the central source (Panaitescu et al. 
2005; Zhang et al. 2005). Another possibility would be to delay the rise of the forward shock 
contribution as a result of viewing angle effects (Eichler and Granot, 2005). These two pro- 
posals would however strengthen the constraint on the efficiency of the prompt gamma-ray 
emission which is already a potential problem for the internal shock scenario (see however 
Fan & Piran, 2006 and Zhang et al. 2006). For the initial steep decay, curvature effects of 
the emitting shell have been invoked (Nousek et al. 2005) while flares are usually explained 
by a late activity of the central source (Zhang et al. 2005; Fan & Wei, 2005). 

In this paper we do not consider the origin of flares but rather focus on the evolution of 
the early X-ray afterglow. We propose that it could be accounted for by a contribution from 
the reverse shock. We develop a simple model which allows us to follow the internal, reverse 
and forward shocks in a consistent way. We compute the energy dissipated in the reverse 
shock and show that, for some specific initial distribution of the Lorentz factor in the flow, 
it is possible to reproduce the succession of the three slopes a±, a 2 and 0:3. We then discuss 
under which conditions part of this dissipated energy can be radiated in the X-ray range, 
providing an alternative explanation for the early X-ray afterglows of GRBs. We also obtain 
the optical emission of the reverse shock and show that chromatic breaks can be observed 
in some cases. 

The paper is organized as follows: in Sect. 2 we present the simplified model we use to 
follow the dynamics of internal shocks. We explain in Sect. 3 how it is extended to include 
the interaction with the environment and we compute the power dissipated in the reverse 
shock. We consider in Sect. 4 the possibility for this power to be partially radiated in the 
X-ray range. We discuss in Sect. 5 the relative importance of the reverse and forward shock 
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contributions and present X-ray and optical afterglow light curves produced by the reverse 
shock alone. Sect. 6 is our conclusion. 



2 THE ORIGIN OF GRB PULSES 

In the context of the internal shock model for the prompt emission of GRBs the pulses 
observed in the light curve are produced when fast moving material catches up with slower 
one previously ejected by the central source (Rees & Meszaros, 1994). This process has often 
been represented by the collision of two "shells" of negligible thickness. However the central 
source probably does not release individual shells but a continuous relativistic outflow with 
a varying Lorentz factor. For this reason the shape of the pulses is largely dominated by hy- 
drodynamical effects (Daigne & Mochkovitch, 2003) while high latitude emission (curvature 
effect) only becomes important at late times. Soderberg and Fenimore (2001) have for exam- 
ple found that the decay of pulses differs from what would be expected if it was controlled 
by the curvature effect alone. A hydrodynamical study of the relativistic flow emerging from 
the central engine therefore appears necessary for a detailed description of the physics of 
pulses but it is naturally quite expensive in computing time (Daigne & Mochkovitch, 2000). 
Fortunately it can often be replaced by a simplified model where the flow is represented by 
a large number of regularly ejected shells which interact by direct collision only (Daigne & 
Mochkovitch, 1998). This neglects pressure waves but this is a good approximation since 
kinetic energy strongly dominates over internal energy of the flow. This approach implies 
to use many shells (from 10 3 to 10 4 ) to represent accurately the distribution of mass and 
Lorentz factor. It is different from the even more simplified description where the numbers 
of shells essentially corresponds to the number of pulses to be produced and where the tem- 
poral profiles are then entirely fixed by the curvature effect (Kobayashi, Piran & Sari, 1997). 

The position Ri of each shell of mass and Lorentz factor Tj is followed as a function of 
time t (in the source frame). When shell % catches up with shell i + la shock occurs at time 
t s and radius R s . The two shells merge and the resulting Lorentz factor after the collision is 
given by 




(1) 



If the the released energy 



4 Genet et al. 

E = rriiTiC 2 + m i+1 T i+1 c 2 - (rrii + m i+1 )T T c 2 (2) 
can be efficiently radiated it will be received by the observer at a time 

^obs = (3) 
C 

and for a typical duration 

AU, = ^ (4) 

under the condition that the radiative time is much smaller than the dynamical time (fast 
cooling regime). The burst bolometric luminosity can then be obtained from the sum of all 
the elementary shock contributions, the dynamical evolution being terminated when all the 
shells have their Lorentz factor decreasing downstream so that no new internal shock can 
form. For an accurate description of the pulse profile at late times, the contribution £(t) of 
each elementary shock must include the curvature effect of the emitting shell which yields 
IE 

^ = At u (I + t^is ( 5 ) 

^obsll -I- Atobs ) 

for t bs < t < t ohs + (1 — cosA#) R s /c, where A6 is the opening angle of the jet here supposed 
to be seen on axis (Granot, Piran & Sari, 1999; Woods & Loeb, 1999). The luminosity in 
a given energy band depends on some additional (and uncertain) assumptions on the post- 
shock magnetic field and Lorentz factor of the electrons which are discussed in Sect. 4 and 

5 while in Sect. 3 we restrict ourselves to the bolometric emission only. 

To produce a single pulse burst (or a pulse as a building block of a more complex 
burst) we have often used in previous works (Daigne & Mochkovitch, 1998, 2000) an initial 
distribution of the Lorentz factor of the form 

t-i/j\ Tmax "I - r m ; n r max r m j n / t \ , . 

T{t) = 2 2 C ° S (^) (6) 

if t < 0.2 tw an d r(t) = r max if t > 0.2 tw! Tmax and r min are the maximum and minimum 
values of the Lorentz factor and tw the duration of the relativistic wind emission (the first 
shell is then ejected at t — and the last one at t — tw)- This Lorentz factor distribution 
where a rapid part of the flow is decelerated by a slower part placed ahead of it, has been 
represented in Fig. la for r max = 200, r min = 50 and % = 10 s. The resulting bolometric 
profile from dissipation by internal shocks is shown in Fig. lb for a total (isotropic) radiated 
energy E rad = 10 53 erg. 

The decline of the pulse after maximum is first controlled by the dynamics of internal 
shocks. This would lead to an asymptotic behavior L(t) oc t~ 3 / 2 if it was not interrupted 
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Figure 1. A single pulse burst, (a): Initial distribution of the Lorentz factor in the relativistic flow as a function of ejection 
time and distance D to the source (in light-seconds) according to Eq.(6) with r max = 200, r m ; n = 50 and % = 10 s (thick full 
line) and to Eq.(7) with <5 = 1 and Tj = 2 (thin full line); (b): Bolometric profile for the distribution given by Eq. (6) (full line) 
together with the corresponding temporal slope a (dashed line). After maximum, the profile is first controlled by the dynamics 
of internal shocks before the curvature effect eventually dominates after t ^ s ~ 20 s. 

at a time r ~ t w when all the fast material of the ejecta has been shocked. Daigne and 
Mochkovitch (2003) have shown that the gamma-ray profiles which can be obtained from 
this first part of the bolometric light curve are generally in good agreement with the early 
decline following maximum count rate in observed GRBs (Ryde & Svensson, 2000). 

After all the ejecta has been shocked the pulse evolution becomes fixed by geometry, the 
contribution of each shocked shell being given by Eq.(5). At large times t 3> r, all the £(t) 
and therefore the global profile L(t) asymptotically behave as t~ 3 . However at early times 
t > r, a steeper decline can be obtained (Nousek et al. 2005) as illustrated in Fig. lb where 
the temporal slope a has been plotted together with the profile. It has a maximum of 3.65 
just at the end of the internal shock phase before relaxing to 3 after a few r. 



3 THE REVERSE SHOCK 
3.1 Physical description 

The profile calculated above corresponds to a "naked GRB" (Kumar & Panaitescu, 2000) 
and would be the only component observed in the absence of external medium. The burst 
environment will however interact with the ejecta, leading to a forward shock propagating 
through the circumstellar medium and a reverse shock sweeping back into the ejecta. We 
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compute below (Sect. 3. 2) the power dissipated in the reverse shock and discuss its possible 
contribution to the early X-ray emission of GRBs. With the initial distribution of the Lorentz 
factor given by Eq.(6) the reverse shock crosses the ejecta in a short time and cannot explain 
an emission lasting for several days. The situation is however very different if a slight change 
is made in the initial distribution of the Lorentz factor. We expect that the central source 
will not stop ejecting relativistic material abruptly at t — t w . We instead propose that T 
will progressively decrease until it reaches a small value (possibly close to unity) at %. Since 
F(t) is given by the ratio E/M of the energy to mass injection rates, a small V can be a 
consequence of (i) a decrease of E, less and less energy becoming available from the source 
to accelerate a given baryon load or/and (ii) a catastrophic increase of M. Case (i) appears 
more natural during the late stages of source activity and has been adopted in presenting 
our results in Sect. 3. 2. 

We have then introduced a new distribution of the Lorentz factor where, for t > 0.5 tw, 
T(t) decreases to a final value Tf 

5 

(7) 

while for t < 0.5%, T(t) is still given by Eq.(6). This modified Lorentz factor is plotted 
in Fig.la for r max = 200, T f = 2, t w = 10 s and 5 = 1. With this new distribution, the 
duration of source activity remains unchanged but the reverse shock will be present for a 
much longer time, until all the ejecta has been decelerated to T ~ Tf . The forward shock also 
remains feeded in energy as slow material from the ejecta is continuously catching up but 
the resulting effect is too small in this case to account for the shallow part of the light curve 
(we assumed that equal amounts of kinetic energy are injected before and after t = 0.5%)- 
To compute the energy dissipated in the reverse shock we had to implement in our shell 
model the interaction with the burst environment. This was done by considering the contact 
discontinuity which separates the ejecta and the shocked external medium. In our simple 
description it is represented by two shells moving at the same Lorentz factor T. The first 
one corresponds to the mass M ej of the ejecta already crossed by the reverse shock, which 
carries a total energy TM^c 2 , and the second to the shocked external medium of mass M ex . 
If the forward shock moves quasi-adiabatically (slow cooling regime), this shell keeps its 
internal energy (since pdV work is neglected in our simple model) so that its total energy is 
TTiM ex c 2 where (r, — l)c 2 is the dissipated energy (per unit mass) in the fluid rest frame. 
Two processes will affect this two shell structure at the contact discontinuity : it will collide 



r(t) = r f + (r, 



r f ) 



i - t/t w 
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either with shells of the external medium at rest, or with rapid shells of the relativistic ejecta 
catching up. This represents both the forward and reverse shock in our simplified picture. 

Forward shock: the interaction with the external medium is discretized by assuming that a 
collision occurs each time the contact discontinuity has travelled from a radius R to a radius 
R' so that the swept-up mass is 

rR' m 

m cx = / 4nr 2 p(r)dr = q— (8) 
Jr 1 

where M = M e j + M ex , p(r) is the density of the external medium (for which we adopted 
either a constant or a stellar wind distribution) and q <C 1 (we take in practice q = 1CT 2 ). 
Writing the conservation of energy-momentum for this collision, we obtain the new Lorentz 
factor r r at the contact discontinuity 



(M ej + M ex r,)r 2 + m ex ri 1/2 



(9) 



(M ej + M ex r 4 ) + 2m cx r 
and also the new Lorentz factor 1^ for internal motions after the collision 

r , = (M ej + M ex r,)r + m cx - M cj r r 

(M cx + m cx )r r • 1 ' 

It should be noted that the above equations assume that material in the burst environment is 

at rest. This neglects the pair- loading process resulting from the initial flash of gamma-rays 

which pre-accelerates the circumstellar medium (Madau & Thompson, 2000; Thompson & 

Madau, 2000; Beloborodov, 2002) out to a radius 

R &cc ~7l^E\% cm (11) 

where -E 7i 53 is the isotropic gamma-ray energy of the flash in units of 10 53 erg. For this 
reason, the deceleration by the external medium will be delayed by 

Atacc~^ = 12< 5 2 3 r 2 - 2 s (12) 

where Y% is the average Lorentz factor of the ejecta in units of 10 2 . Therefore the initial 
dynamical evolution will be that of a naked GRB but the effect will last more than one 
minute only for the most extreme bursts with -E 7i 53 > 10. 

Reverse shock: as the Lorentz factor at the contact discontinuity decreases, new shells from 
the ejecta become able to catch up. Writing again the conservation of energy-momentum for 
these collisions, we obtain the change in Lorentz factor 



(M ej + Me X r i )r + m C j7 cj n ] : 
(M ej +M cx r i ) 7ej +m cj r 



(13) 
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Figure 2. Dissipated power as a function of observer time during reverse shock propagation when the Lorentz factor is given 
by Eq.(7) with T f = 2 and 5 = 1. (a): wind case with A, = 0.5 (dashed line), A* = 0.1 (full line) and A* = 0.05 (dotted line); 
(b): uniform density case with n = 1000 (dashed line), 10 (full line) and 0.1 cm -3 (dotted line). In each panel the thin full line 
represents the naked burst. 

and the related dissipated energy 

£diss = (M cj + M^T^Tc 2 + m cj7cj c 2 - (M cj + m ej + M^r^iV (14) 

where m e j and 7 C j are respectively the mass and Lorentz factor of the colliding material from 
the eject a. 

3.2 The dissipated power 

Using this simplified model for the interaction of the ejecta with its environment we can 
describe the deceleration of the front shell and the propagation of the reverse shock. We 
have obtained the dissipated power in the reverse shock for different burst environments 
(uniform medium or wind). In the wind case, we considered three values of the parameter 
A*: 0.5, 0.1 and 0.05 (such as p(r) = 5 10 11 A*/r 2 g.cm -3 with A* = 1 for a wind mass 
loss rate M w = 10~ 5 M .yr _1 and a terminal velocity = 1000 km.s -1 ). In the constant 
density case, we also tried three values of n: 1000, 10 and 0.1 cm~ 3 . The resulting profiles are 
shown in Fig.2 for T f = 2 and 5 = 1 in Eq.(7) but we checked that they remain essentially 
unchanged when Tf is varied between 1 and 10 and S between 0.5 and 2. 

The curves in Fig.2 show a striking resemblance with the early X-ray afterglows observed 
by SWIFT. After about 100 s the reverse shock component dominates over the tail of the of 
the burst prompt emission computed in the last section. At late times the decline follows a 
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constant slope a ~ 1.5. The shape of the intermediate region is most sensitive to the density 
of the burst environment. At high density it is nearly suppressed, the constant slope a ~ 1.5 
following directly the initial steep decrease. Conversely at low density, it can become com- 
pletely flat and even fall to a temporary minimum. For comparison, we have also represented 
in Fig. 2 the profile corresponding to the naked GRB. Without an external medium and for 
the distribution of Lorentz factor given by Eq.(7) (with T f = 2 and 5 = 1) we find that a 
large fraction of the ejecta (all slow material with T < 140) remains unaffected by internal 
shocks. With an external medium the reverse shock propagates throught this material which 
produces the additional power at late times. 



To elucidate the behavior of the reverse shock contribution we have considered the fol- 
lowing simplified case which can be handled analytically: the ejecta is supposed to be made 
of a rapid single shell of mass Mq and initial Lorentz factor To (representing the fast material 
where the prompt emission takes place) followed by a slower tail of the form 

r(M) = r f + (r - r f ) (15) 

where T f is the final Lorentz factor at the end of the tail and M s is the total mass of the 
slow material (M = corresponds to the last emitted shell and 5 allows to vary the tail 
shape). Notice that this expression of T(M) directly results from Eq.(7) if M is constant. 
Such a distribution of T skips the prompt phase and the resulting dissipated power P diss 
comes from the reverse shock only. 

The reverse shock contribution is maximum at a time close to the deceleration time of 
the front shell 



W — Q d p C 2 with R dcc 

ZC 1 n 



M (3 



(16) 



A-kATq 

where we have used the notation p = Ar~ s with A = p and s = for a uniform medium 
and s = 2 for a stellar wind. A full analytical solution can be obtained for the reverse 
shock contribution but we only derive below its asymptotic behavior assuming that the 
front shell essentially follows the Blandford-McKee solution, i.e. it is only weakly affected 
by the additional energy coming from the slow material progressively catching up. We have 
checked this approximation with the numerical simulations and it is satisfied to an accuracy 
of about 25%. 

We can write the power dissipated in the reverse shock as 
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dMdT^ 2 

p «" = -WTi Tec (17) 

where t is the observer time when a shell of Lorentz factor T catches up with the front shell. 
The fraction e of the incoming material kinetic energy dissipated in the collision can be 
obtained from Eq.(13) and (14) with m C j <C M ej + M ex Fj, leading to 



1 - 



r 



(18) 



(20) 



where Tf s is the Lorentz factor of the front shell given by the Blandford-McKee solution 

r fs ^r (-^V A (19) 

with A = = 3/2 (resp. 1/2) for a uniform medium (resp. a stellar wind). Using dR/dt = 
2c r| we then get the relation between shock radius and observer time 
t _1 / R \ 2A+1 

^dcc 2 A + 1 V-Rdcc/ 

With our assumed distribution (Eq.(15)) of the Lorentz factor in the slow material which is 
steadily increasing outwards, each shell moves independently at the constant Lorentz factor 
T until it catches up with the forward shock. We moreover neglect the fact that the slow 
material is emitted over a certain duration and we write the position of each shell as a 
function of observer time as 

R ' < T ^ . (21) 



Rdec tdcc VTo 

Notice that if the Lorentz factor is not initially monotonic in the slow material, internal 
shocks will take place which, when they are completed, will leave a new distribution of T, 
now steadily increasing outwards. For observing times long compared to the time of internal 
shocks Eq.(21) will therefore still hold. 

Eliminating the radius between Eq.(20) and (21) gives the time when a shell of Lorentz 
factor T catches up with the forward shock 

2A+1 

1 (2A + iWry * _ (22) 



^dcc v r 

Now from Eq.(19), (20) and (22) we can obtain the Lorentz factor Tf s of the forward shock 
when the slow shell of Lorentz factor Y catches up 

r fs = r (2A + 1)- 1 / 2 . (23) 

To end with a simple power law expression for the dissipated power we write 



<ir st„ Vr„ 
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Figure 3. Synthetic X-ray light curves in the 0.3 - 10 kcV range for the models of Fig. 2. The adopted post-shock energy 
redistribution parameters are e e = eg = 1/3 and £ = 10~ 2 . The assumed burst redshift is z = 2; The different cases considered 
in Fig.2 are represented by the same full, dotted and dashed lines. 



i.e. we use Eq.(15) with Tf = which is obviously uncorrect but does not change the behavior 
of the solution in the relativistic regime. Using Eq.(22), (23) and (24) we finally get 
r M s c 2 / t y 



P diss (t) = $(A,5) 



where 



& tdec V ^dec / 



S(A,*) = ^ [l-(2A + l)- 1/2 ; 2 



l-^(l+4A) 
X (2A + 1) 2«(2A+1) 



(25) 



(26) 



For 5=1 and the two values of interest for A, Eq.(25) becomes 

P / i \ -7/4 

- = 6.610- 2 ( — ) (A = 3/2) 

r M s c 2 /t dec \t dc J 

t \- 3 / 2 



1.510- 2 ( — ) (A = 1/2) 

^ ''dec ' 



(27) 



which is in good agreement with the asymptotic behavior of the light curves shown in Fig.2. 



4 CAN THE REVERSE SHOCK CONTRIBUTE IN X-RAYS? 

Despite their similarity with the SWIFT observations, it must remain clear that the profiles 
shown in Fig.2 only trace the power dissipated in the reverse shock. With the assumptions 
ordinary made to compute the reverse shock contribution in GRBs it should manifest itself 
mainly in the visible/IR range (Sari & Piran, 1999). Moreover most of the emission would 



12 Genet et al. 

generally take place in the slow cooling regime so that the observed light curve will not 
necessarily trace the instantaneous energy release. 

We therefore investigated whether, under some specific conditions, a substantial fraction 
of the dissipated power can be (i) radiated in the X-ray range and (ii) in the fast cooling 
regime. If the reverse shock contribution originates from synchrotron radiation of shock 
accelerated electrons, the characteristic synchrotron energy and cooling time behave as 

E s oc BT 2 e and t s oc B~ 2 T~ l (28) 

in the rest frame of the shocked material. Both the post-shock magnetic field B and typical 
electron Lorentz factor T e have therefore to be large to produce an emission at high energy 
and on a short time scale t s < t& yIl . An estimate of T e is usually obtained assuming that a 
fraction e e of the dissipated energy is injected into a fraction ( of the electrons so that 

r e ^^e (29) 

where m p and m e are the proton and electron masses and ec 2 is the energy dissipated per 
unit mass in the comoving frame. Similarly the post-shock magnetic field can be expressed 
as 

B = (8ne B pec 2 ) 1/2 (30) 

where p is the comoving density and e# the fraction of the dissipated energy tranferred 
to the magnetic field. To have large B and T e values behind the shock we first supposed 
that a complete equipartion is established between the electronic, magnetic and baryonic 
components so that e e = €b = ebaryon = 1/3. We also assumed that only a small fraction 
C ^ 1CT 2 of the electron population is accelerated in the shock. Adopting ( = 1CT 2 increases 
T e by a factor of 100 and hence E s by a factor 10 4 and decreases t s by 10 2 compared to the 
standard ( = 1 case. 

The possibility to have only a small fraction of electrons being accelerated has already 
been considered by Bykov & Meszaros (1996) and also by Eichler & Waxman (2005) in the 
context of GRB afterglows. They showed that ( is not well constrained by the observations 
and, even if ( ~ 1 appears slightly favored, they included the whole interval m e /m p < ( < 1 
in their analysis. In internal shocks, which are very similar to the reverse shock (both take 
place in the burst ejecta and are mildly relativistic) a large e e is required to maintain a 
reasonable global efficiency since the fraction of the total energy dissipated by internal 
shocks hardly exceeds 10%. A small ( is also favored to insure that the emission takes place 
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in the gamma-ray range as shown by Daigne and Mochkovitch (1998) and more recently by 
Lee et al. (2005) in the context of the short hard burst GRB 050509b. 

Examples of synthetic light curves in the XRT band 0.3 - 10 keV are shown in Fig. 3 
for the cases already considered in Fig. 2. They have been obtained with e e = e# = 1/3, 
( = 10~ 2 , a slope p = 2.5 for the electron energy distribution and an assumed redshift 
z = 2, typical of the SWIFT burst population. Especially in the wind case, they seem able 
to reproduce many of the observed XRT light curves. Conversely in the uniform density case 
we often obtain a depressed minimum followed by a bump rather than a continuous shallow 
evolution 



5 DISCUSSION 

Our proposal to explain the early X-ray afterglow of GRBs by a contribution of the reverse 
shock relies on three well defined assumptions: (i) the Lorentz factor of the material ejected 
at late times by the source has to decrease to small values, T f < 10; (ii) the shock dissipated 
energy must be transferred to only a small fraction of the electron population; and (Hi) the 
forward shock contribution should lie below that of the reverse shock, at least during the 
first hours following burst trigger. 

This last condition requires an ineffective transfer of energy to electrons (e e ^ 10 -2 ) 
or/and magnetic field (eb ^ 10~ 5 ) in the material crossed by the forward shock. Difficulties 
to produce a sufficiently large magnetic field extending over the emitting region of GRB 
afterglows has for example been recently emphasized by Milosavljevic & Nakar (2006). Then, 
if the reverse shock dominates in X-rays, what is the situation in the visible? We have 
checked that in most cases, taking small values of e e or/and e# in the forward shock, equally 
implies that the reverse shock dominates in the visible. The consistency of our proposal must 
therefore be checked not only with X-ray observations but also at lower wavelengths. 

To better understand the multiwavelength behavior of the reverse shock contribution, we 
have computed the peak flux F max and the characteristic synchrotron and cooling frequencies 
v m and v c (Sari, Piran & Narayan, 1998). These three quantities depend on t (observer time), 
N e (total number of shock accelerated electrons), B (magnetic field in shocked material), 
T e (typical electron Lorentz factor) and V (Lorentz factor of the emitting material) in the 
following way 
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F max oc TBN e 



v m <xTT\B (31) 

We consider their temporal evolution in the asymptotic regime already described in Sect. 3. 2, 
assuming a wind environment. The evolution of T is given by Eq.(22) which can be reex- 
pressed as 

/ t \~ 1/4 / t \" 1/4 

r(t) = r ( — ) =ii9r 2 ( — ) (32) 



where T 2 = r /100. Eq.(18), (23) and (29) show that T e reaches a constant value 

r e = 4.31(T 2 ^^ = 79^^ (33) 

where we have added the normalizing factor (p being the slope of the relativistic electron 
distribution) which was not present in Eq.(29). For the magnetic field, instead of Eq.(30) it 
is easier to use the continuity of the energy density at the contact discontinuity which yields 

B = (32ne B c 2 A) 1 / 2 ^ (34) 

where A is the wind constant such as p(R) = A/R 2 (A = 5 lO 11 ^). With Eq.(20) for R and 
Eq.(32) for T we get 

B (t) = 3 10 4 (eBA ; )V2 ( —) G (35) 



^dcc 1^2 V td. 

Finally, the number of accelerated electrons can be obtained from Eq.(24) which, for S — 1 
gives 

t \- 4 / 4 



1 - 1.19 



(36) 



T m p c 2 

where E s = |r M s c 2 is the total energy in the slow material. From Eq.(32), (33), (35) and 
(36) the expressions for F max , u m and u c can be computed 

F max = 1.4 10 8 ^ m^A^l x I mJ 

u m = 9.15 10 16 {eU^' 2 (f f (g) 2 x \ Hz (37) 
v c = 8.2 10 8 jgZfa x t^ 2 Hz 

with E 53 = E s /10 53 erg and where the expression for F max has been written in the limit 
t S> tdcc- Compared to the forward shock case, it can be seen that F max oc t' 1 and that v m 
decays less rapidly (as t^ 1 instead of t~ 3 / 2 ). For a wind environment, the cooling frequency 
has the same power law dependence, v c oc t 1 ^ 2 . From these expressions the flux can be 
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computed for the different possible radiative regimes (Sari, Piran & Narayan, 1998), the 
results being given in Appendix A. 

Let us for example take the following values of the parameters: T 2 = 1, e e = e# = 0.33, 
( = 10~ 2 , p = 2.5, A* = 0.5 and t dcc = 100 s. Then, the transition from fast to slow cooling 
occurs at t — 1.3(1 + z) day. Now adopting 1 keV and 2 eV as typical energies for the X-ray 
and visible bands (i.e. v x = 2.4 10 17 Hz and v v = 4.8 10 14 Hz) and a redshift z — 2, it appears 
that after only a few seconds vx becomes larger than v m and then remains larger than v c 
in the slow cooling regime. The corresponding temporal slope is ax = (2p + l)/4 = 1.5. At 
the visible frequency, we initially have v c < u v < v m and therefore ay = 0.75. The visible 
frequency crosses v m at t — 2.6 h (in the fast cooling regime) and then v c (in slow cooling) 
at very late times. A break from a v = 0.75 to 1.5 is expected at t — 2.6 h. 

Since these predicted slopes are only valid in the asymptotic regime where t 3> tdec we 
have performed a numerical simulation with the burst parameters given above except for 
the fraction ( of accelerated electrons which is varied between 0.003 and 0.03. We assume 
in addition that E 53 = 1 and adopt a rest frame reddening Ay = 0.5 in the burst host 
galaxy. The resulting X-ray and visible light curves are shown in Fig. 4. For ( = 3 10~ 3 and 
10~ 2 they exhibit chromatic breaks. The break in X-rays is a consequence of the dynamics 
of the reverse shock (it is already present in the bolometric light curve) while the break in 
the visible is a spectral break (when vy crosses is m ). The cases with ( = 3 10~ 3 , 10~ 2 and 
3 10~ 2 are very similar to the early afterglow light curves of respectively GRB 050802, GRB 
050922c and GRB 050801 (see Panaitescu et al, 2006 and Panaitescu, 2006). 

The subsequent evolution of the afterglow will depend on the behavior of e e and e# 
in the forward shock. If they increase enough with time the forward shock contribution 
will eventually dominate but the moment of the transition is difficult to estimate in the 
absence of any reliable physical model for the possible variations of the shock microphysics 
parameters. If the forward shock takes over after about one day, the multiwavelength fits of 
GRB afterglows obtained in the pre-SWIFT era will remain valid but the early afterglow 
will be explained by the reverse shock. At the transition, a change of slope or the presence of 
a bump may however be expected. While such accidents have been observed in some bursts 
they do not seem to be a generic feature of GRB afterglows. 

A much more radical point of view can still be adopted: it would to suppose that in some 
cases the forward shock never takes over so that the afterglow is entirely produced by the 
reverse shock! The results shown in Fig. 4 seem to indicate that this possibility should not be 
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times (hours) times (hours) times (hours) 



Figure 4. Early afterglow light curves produced by the reverse shock for e e = eg = 0.33, p = 2.5, T2 = 1, A* = 0.5, t^cc = 100 
s, Ay = 0.5 and, from left to right, £ = 310~ 3 , 10~ 2 and 310~ 2 (see text for details). The full (resp. dotted) line is the 
X-ray (resp. the visible) afterglow. Compare these results to the early afterglow light curves of respectively GRB 050802, GRB 
050922c and GRB 050801 as shown in Panaitescu (2006). 

excluded a priori even if, clearly, considerable work will be needed to confirm it. As for the 
forward shock hypothesis it will have to be confronted to a large amount of multiwavelength 
afterglow data and show that it can provide a consistent picture for their interpretation. 

6 CONCLUSION 

We have developed a simplified model which enabled us to follow simultaneously the dy- 
namics of the internal, external and reverse shocks in GRBs. We were mainly interested by 
dissipation in the reverse shock when the Lorentz factor in the material which is ejected at 
late times by the source decreases to small values, T f < 10. The propagation of the reverse 
shock then extends over quite a long time needed to decelerate the fast moving part of the 
ejecta down to T ~ T f . We have obtained the dissipated power as a function of observed 
time for different burst environments (wind or constant density). Its evolution shows a stri- 
king resemblance with the early afterglow light curves observed by SWIFT, especially in the 
wind case. However the reverse shock contribution is normally expected at low energy and 
to appear in X-rays it requires a transfer of the dissipated power to only a small fraction 
(C ^ 1CT 2 ) of the electron population. If this is possible, SWIFT XRT observations could 
be better explained by the reverse shock than by the standard afterglow produced by the 
forward shock. 

We have also computed the optical emission from the reverse shock. The comparison 
with the X-ray light curve often reveals the presence of chromatic breaks during the first 
hours. Such breaks have been observed and are difficult to explain with the standard model 
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where the afterglow comes from the forward shock. We have finally even proposed that in 
some cases the entire afterglow could be produced by the reverse shock. We fully understand 
that, to be validated, this non standard view still has to show that it can successfully explain 
multiwavelength observations of a reasonable sample of GRB afterglows. We aim to perform 
these necessary tests in a work in preparation. 



APPENDIX A: 

Using Eq.(37) for F max , v m and v c we give the expressions for the flux at a given frequency 
in the fast and slow cooling regimes: 

Fast cooling 
1) v < v c 

F F f^ 1/3 



V, 



^■ '^r^y ^x^ mJ (A1) 



2) V c < V < V m 

V \ -V2 



F = F — 

1 v 1 max 1 



7 7 IP 3 (j±^! (^ 53 ) ^cc -1/2 ,-3/4 T 

3) v > v m 



A v A max I 1 i 



1.3 10 4 x 0.36-/ 2 (1 + *t P/2 ((E 53 ) {e B A^tt c 

IJ- 



l-p/2 
28 



S7ou; cooling 



, V \ V3 

-^max 



210 s 



-^28 r 2 

"2/3 / \ -2/3 



(A2) 



x (f)" 1 ^)" 1 ^ 2 -^ -J (A3) 



x I (0) & x ^ mJ (A4) 
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2) v m < v < v c 

F = F 



max 



l-i 

V \ 2 



1.410 s x 0.36 2 



Dno To 



y 28 1 2 

x (|) P " 1 (^T) P " 1 "- ,/2xr41 mJ (A5) 

3) v > u c 

= 1.3 10* x 0.36'/ 2 il±i|^! (( b 53) (es . 4>) ^ t V4 

x (l)" 1 ^)" 1 "-^'"^ mJ (A6) 

In all these expressions the frequency (in observer frame) is in unit of 10 17 ' 4 Hz, corresponding 
to 1 keV. 
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